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ABSTRACT 

We report the discovery of a systematic constant time lag between the X-ray and radio flares of 
the gamma-ray binary LSI +61 303, persistent over long, multi-year, time scale. Using the data of 
monitoring of the system by RXTE we show that the orbital phase of X-ray flares from the source varies 
from <j)x — 0-35 to <j)x — 0.75 on the superorbital 4.6 yr time scale. Simultaneous radio observations 
show that periodic radio flares always lag the X-ray flare by A<px-R — 0.2. We propose that the 
constant phase lag corresponds to the time of flight of the high-energy particle filled plasma blobs 
from inside the binary to the radio emission region at the distance ~ 10 times the binary separation 
distance. We put forward a hypothesis that the X-ray bursts correspond to the moments of formation 
of plasma blobs inside the binary system. 
Subject headings: binaries - gamma-rays, X-rays, radio: individual: LSI +61 303 



1. INTRODUCTION. 

LSI +61 303 is one of several high- mass X-ray binaries 
which emit hi gh-energy (GeV -TeV) 7-r ays (see e.dAlbert 
et al.| ( |2006D ; |Acciari et al.| ( |2008[ ); |Abdo et al.j ( |2U09p j. 
Although the nature of high-energy activity in high-mass 
X-ray binaries is not clearly understood, a common hy- 
pothesis is that it is related to existence of relativistic 
particle outflow. 

LSI +61 303 is known to be variable on different time 
scales. The orbital period is P or & — 26.496 d ( |Gre- 



2002|. Zero orbital phase </> — corresponds 
2 443 366.775 + nP or t), periastron occurs at 



orb 



gory 
to io'_ 

(/> = 0.275 (Ara gona et al . 2009). The superorbital pe 
riod is P so =1667d (Gregory 2002} . The zero superor- 
bital phase $ = is T , so = 2 443 366.775 + kP so . In 
LSI +61 303 the high-energy particle outflow is directly 
observed in the radio band, where angular resolution is 
sufficient to resolve the source and detect variations of its 



rather a variable morphology outflow filling a region the 
size ~ 10 2 — 10 3 times larger t han the binary separation 
distance (D hawan et al. ||2006[) . 

The radio signal could not be used to trace the out- 
flow down to the production site inside the binary orbit, 
because of the free-free absorption in the den se stellar 
wind environment (e.g. ( Zdziarski et al.|[2010| )). To un- 
derstand the nature of the high-energy particle carrying 
outflow one has to use complementary high-energy data 
in X-ray and/or 7-ray bands. 

The X-ray and 7-ray emission from the sys tem is 
known to be variable on the orbital time scale ([H arri- 
son et al.||2QQQ| |Chernyakova et al.||2QQ6J |Zhang et aT 



morphology on the orbital period time scale (Paredes et 
al.|19 98; Mas si et al.|2004| ). The observed morphological 
changes indicate that the outflow is, most probably, not 
a jet with a well defined position angle on the sky, but is 



2010p . A"study of the orbital modulation of the X-ray 
and 7-ray signal was missing up to recently due to the 
absence of systematic monitoring of the source on many 
orbit (year) time scale. Such a monitoring has recently 
became possible in the GeV 7-ray band with the start of 
operation of the Large Area Telesco pe (LAT) on board 
of Fermi satellite ( Abdo et al.||2009 ) which has now col- 
lected ~ 3 yr of data. In the X-ray band, a dedicated 
multi-year monitoring campaign was done with RXTE 
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( [Smith et al.||2009||Torres et al.|2010||Li et al.|2011| >. As 
a result, both the average periodic modulation pattern 
and orbit-to-orbit variations of X-ray and 7-ray emission 
from the source were established. 

In the X-ray band the source is known to exhibit one 
flare per orbit, on average preceeding the radio one. The 
GeV band lightcurve during the first year o f LAT obser- 
vations also exhibited a one-flare pattern (Abdo et al. 
2009). Both X-ray and GeV band light curves exhibit 
large orbit-to-orbit variations so that the systematic pe- 
riodic variability is often washed out by an erratic be- 
haviour of the source. The origin of the X-ray and radio 
flares as well as the relation between the flaring activity 
of the source in different energy bands is not well under- 
stood. 

The orbital phase of the periodic flares drifts on a 
super-orbital time scale by half-an-orbit from (j) r ~ 0.5 to 
(j) r ~ 1 (Gregory 2002). Such a drift is difficult to explain 
in scenarios based on various types of precessio n, where 
one e xpects a drift by a full orbit < <p r < 1 (Gregory 
2002). An alternative possibility for the explanation of 
the 4.6 yr time scale is the build up and decay of the 
equa torial disk around the massive Be star in the system 
flZamanov et al.||1999| ). 

A new insight in the nature of the 4.6 yr periodic- 
ity/variability might be given by the study of the changes 
in the behaviour of the X-ray and 7-ray emission on this 
time scale. Here we report such a study based on the 
analysis of the monitoring of the source with RXTE, IN- 
TEGRAL and Fermi. The X-ray/7-ray data are comple- 
mented by the contemporaneous radio monitoring data. 

2. DATA ANALYSIS. 

In our analysis we consider radio, X-ray and 7-ray data 
on LSI +61 303 collected over the period from 2007 to 
2011. 

LSI+61 303 was regularly observed by the PC A (Bradt 



et al. 1 1 1993 ) instrument on board of RXTE (one 
2 — 3 ks observation every 3-5 days) since Sept. 2007 
up to Aug. 2011. During this period the set of operat- 
ing detectors changed from observation to observation. 
Only detector number 2 (PCU2) operated all the time. 
For uniformity we consider only the data from PCU2. 
To measure the source flux for an individual observation 
we extracted an appropriate mean spectrum using the 
Standard-2 mode data. For background subtraction we 
used the latest PCA backgrou nd model for faint sources 
(see e.g. Jahoda et al. (2006) and PCA web-page). For 
the wide energy band (c>-20 keV) the backround estima- 
tion for pca2 is better than 0.1 mCrab. There are no 
known instrument related timescales close to orbital and 
superorbital periods. Most of the data analysis was per- 
formed using the HEASOFT 6.11 software package. 

LSI +61 303 was several times in the field of view of IN- 
TEGRAL imager ISGRI during the Galactic plane scans 
and pointed observations. In our analysis we use all pub- 
licly available data for the period from the January 2003 
until October 2010. We process the data using the I NTE- 
GRAL Offline Scie nce Analysis (OSA version 9.0) flCour- 
voisier et aL|2003| ). 

In the 7-ray band we use the data of Large Area Tele- 
scope (LA T) on board the FERMI satellite flAtwood et 
al. ||2 009) collected over the period from August 2008 to 
October 2011. We process the data using Fermi Science 







■ / / 

■ flares 


^ 


PCA/ 


RXTE 






• — r 



— L. 






^_ 


^r 


-N 


^-Li>^ 


\ : 






-MM5H3E0HB 


^ 






, , , i , 


. . i , , , 


■ ■ 





Superorbital phase 



Fig. 1. — 3 - 20 keV flux from LSI +61 303 as a function of 
the orbital vs. superobital phase. The color scale is expressed in 
mCrab units. 

Tools version 9r23pl [|] The likelihood analysis is per- 
formed in a 10° radius region around the source taking 
into account all sou rces from the two-year LAT catalog 
(|Abdo et al.||2011|. Spectra of individual sources are 
rlav 



modeled as power laws with the slope fixed to the best 
fit value from the two-year catalogue and free normaliza- 
tion. 

The radio flux density was measured at 15.4 GHz with 
the Ryle Telescope and its successor, the AMI Large Ar- 
ray (Cambridge UK). The Ryle Telescope used a band- 
width of 750 MHz, and the AMI uses 4 GHz. Pooley & 
Fender (1997) describe the operation of the Ryle Tele- 
scope in observations of this type; a very similar tech- 
nique is used for the AMI. Observations are typically 30 
- 60 min in duration. 

3. RESULTS 

Fig. [l] shows the 3-20 keV X-ray flux from the source 
as a function of the orbital and superorbital phase (as 
seen by RXTE/PCA). From Fig. [l] one could see that 
the source exhibits on average one episode of increased 
X-ray activity per orbit. A regular increase of X-ray ac- 
tivity in the phase interval A<frx — 0.2 — 0.4 is accompa- 
nied by random variations of the source flux, with short 
flares appearing on the time scales Tfi ares ^x <C A(j)xP rb- 
The average source flux in 3-20 keV energy range during 
the active/quiet part of the orbit indicated by the white 
dashed lines on Fig. fl] is about 1 mCrab/0.5 mCrab. 
Typical error of the flux measurment, ~ 0.1 mCrab, is 
dominated by the uncertainity of the PCA background. 
The figure also has the color scale expressed in mCrab 
units. 

The average phase of the X-ray activity period (\>x 
varies on the superorbital time scale. From Fig. [I] we 
find that <\>x exhibits a systematic drift from <fix — 0.35 
to (j>x — 0.75 within one superorbital cycle. Such a drift 
is similar to the systematic drift of the phase o f the pe- 
riodi c radio flares from <pR ~ 0.5 to <pR c± 1 (Gregory 
[2002] ). 

Evolution of the orbital variability of the source in hard 
X-rays (20 - 60 keV) on the superorbital time scale, ob- 
served by INTEGRAL, is shown in Fig. [2] One could see 

1 http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 
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Fig. 2. — Averaged orbital variability of the hard X-ray flux (20 
- 60 keV) from LSI +61 303 for the < $ < 0.5 (black solid 
crosses) and 0.5 < <E> < 1 (red dotted crosses) superobital phases. 
The corresponding exposures (from left to right) are 525, 664 and 
452 ksec for the < $ < 0.5 and 298, 132 and 207 ksec for the 
0.5 < $ < 1. 
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Fig. 3. — Radio flux from LSI +61 303 as a function of the orbital 
and superorbital phases. The color scale is expressed in mjy units. 

that, similarly to the 3-20 keV range, the maximum flux 
happens during the orbital phase 0.25 < <\> < 0.5 in the 
superorbital cycle phase 0.5 < $ < 1. The maximum 
becomes wider and shifts toward 0.25 < <px < 0.75 in 
the superorbital phase < $ < 0.5. 

Observations of the systematic drift of the phase of the 
radio flares from the source reported by Gregory (2002) 
were performed several superorbital cycles before the X- 
ray monitoring campaign by RXTE. To verify the long- 
term stability of the range of the shifts of (J)r over many 
superorbital cycles we use the data of monitoring of the 
source in the radio band which are contemporaneous with 
the RXTE monitoring campaign. Fig. [3] shows the ra- 
dio flux of the source as a function of the orbital and 
superorbital phases. Colorbar shows flux measured in 
mjy. Co mparing Fig. [3] with the equivalent figure from 
Gregory| ( |2002[ ), we find that the overall drift pattern of 
the phase ol the radio flare remained stable over sev- 
eral superorbital cycles. The same drift from the phase 
<j) R ~ 0.55 to (j)R c± 0.95 is observed also in the radio data 
contemporaneous with the RXTE monitoring campaign. 

Comparison of X-ray and radio superorbital variability 
patterns is shown in Figure [3J The average phase of the 
X-ray activity period always preceeds the phase of the 
radio flare by A(j)x-R — 0.2, which corresponds to the 
time delay AT X -r = A(J) X -rP ~ 5.3 d. 

Contrary to the X-ray and radio bands, the super- 



FlG. 4. Very high energy (E>100 MeV) flux from LSI +61 303 
as a function of the orbital and superorbital phases. The color 
scale is expressed in mCrab units. 

orbital modulation pattern is not clearly visible in the 
7-ray band. Fig. [4] shows the source flux in the 0.1- 
10 GeV energy band plotted as a function of the orbital 
and superorbital phase, similarly to Figs. [I] and [3] Col- 
orbar shows flux measured in mCrabs, typical error is 
about 10% of the flux. Long-term source behaviour of 
the source in the GeV band is puzzling. Orbital modu- 
lation was clearly observable at the beginning of Fermi 
observations at the superorbital phase 6.8 < <£> < 6.9. 
The phase of the maximum orbital modulation of the 
GeV flux in this superorbital phase range was close to 
the phase of the X-ray activity. However, in the time pe- 
riod following the superorbital phase <£ ~ 6.9 a clear or- 
bital modulation pattern disappeared (see also Hadasch 
[1 ( |2011| )). Further study of the source on the time scale 
of several superorbital cycles is needed to clarify the re- 
peatability of the observed appearance / disappearance 
of the orbital modulation pattern and its relation to the 
overall 4.6 yr activity cycle of the source. 

4. DISCUSSION. 

A constant time delay between the drifting orbital 
phases of X-ray and radio flares could be naturally ex- 
plained if one takes into account that radio and X-ray 
emission originate from different regions. The radio emis- 
sion is prod uced at large distance from the binary, Dr > 
5 x 10 13 cm ( [Zdziarski et al.|2010| ) while the X-ray flux is 
most probably produced at shorter distances of the order 
of the binary separation 3 x 10 12 cm < Dx < 10 13 cm. 
Assuming that injection of high-energy electrons respon- 
sible for the X-ray and radio flares happens in the same 
event in the binary, one could attribute the time delay 
between the X-ray and radio flares to the time-of-flight of 
the high-energy particle filled plasma to the radio emis- 
sion region. The cooling times of high energy particles 
are long enough to allow them to travel to the radio emis- 
sion region. This gives an estimate of the plasma out- 
flow velocity vr ~ Dr/ATx-r — 10 8 cm/s, which is 
in good agreement with the asymptotic velocity of the 
stellar wind from t he massive star in bot h the polar and 
equatorial regions ( Zdziarski et al.|2010 ). 

The constant time delay between the X-ray and radio 
flare phases suggests the following scenario of production 
of the periodic radio flares. Once per orbit, an event of 
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interaction of the compact object with the stellar wind 
leads to injection of high-energy particles into the stel- 
lar wind. The high-energy particles mixed into the stellar 
wind escape from the binary system with the stellar wind 
velocity. A radio flare occurs at the moment when the 
portion of the stellar wind filled with high-energy elec- 
trons reaches the distances D ~ Dr at which the system 
becomes transparent to the radio waves. 

High-energy electrons, responsible for the radio- 
to- X-ray emission, are held in the plasma out- 
flow by the magnetic field. Electrons producing 
synchrotron emission in the radio band have en- 
ergies E e ~ 10 [5/1 G\~ 1,2 [v R /lQ GHz] 1/2 MeV, 
where vr is the frequency of the radio synchrotron 
emission. The main cooling mechanisms for such 
electrons are synchrotron and/or inverse Compton 
emission with the characteristic cooling time scales 



t s - 0.7 [5/1 G\~ 6/z [u R /10GEz]~ L/z yr and t IC ~ 

6.5 [D/10 13 cm] d. Unless magnetic field in the sys- 
tem is much higher than B ~ 10 G, the synchrotron and 
inverse Compton cooling time scales are not shorter than 
the time- of- flight from inside the binary orbit to the ra- 
dio emission region. Thus energy losses do not prevent 
electrons injected in the binary system from traveling to 
the radio emission region on a time scale of several days. 
The X-ray emission is most probably produced via 
the synchrotron mechanism by electrons of the energies 

E e ~ 100 [5/10 G]" 1/2 [E x l$A keV] 1/2 GeV. This con- 
jecture is supported by the observations of fast variability 
of X-ray emission on the time scales t ~ 10 s (Smith et 



al. 2009) which is comparable to the synchrotron cool- 



mg time of 100 GeV electrons in the B ~ 10 G magnetic 
field. The absence of any obvious break / cut-off features 
in the keV-GeV source spectrum with a maximum in the 
~ 10 — 100 MeV range is in favour of interpretation of 
the entire keV-GeV bump as a single spectral compo- 
nent. In such a model the GeV band emission is due to 
the synchrotron emission by the 10-100 TeV electrons. 

The phase of the 7-ray flare is close to the phase of 
the X-ray flare (Fig. fl] and pi), at least during a part 
of the superorbital cycle 6.8 < <£ < 6.9. It is natural 
to identify the phase of the X-ray/7-ray flare with the 
moment of formation of high-energy particle outflow in- 
side the binary. Long term evolution of this phase could 
not be followed in the 7-ray band, because of the large 
width of the activity period after the superorbital phase 
<£ > 6.9. In the X-ray band the width remains finite 
all over the superorbital cycle. Difference in the su- 
perorbital modulation pattern in X-rays and 7-rays is, 
most probably, related to different regimes of accelera- 
tion/propagation/cooling of 10 MeV and 10 TeV elec- 
trons. 

The phase of the X-ray activity is close to the phase of 
the periastron of the binary orbit (j) per ~ 0.3 at the super- 
orbital phase $ c± 0.5 (Fig. IT]) when the duration of the 
activity period is shortest. The phase of activity gradu- 
ally shifts to the post-periastron interval (fix > 0.3 over 
the superorbital period and almost reaches the phase of 
the apastron <fi ap ^0.8 toward the end of the superor- 
bital cycle so that the X-ray emission is always delayed 
with respect to the apastron. Shift of fix is accompanied 
by the increase of the width of activity period. 



A possible explanation for such behaviour could be 
found in a scenario in which the 4.6 yr superorbital cycle 
is interpreted as the cycle of gradual buildup and decay 
of the equatorial disk of the Be star. At $ ~ 0.5 the 
equatorial disk is weak. The phase of the closest en- 
counter between the disk and the compact object is the 
phase of the periastron. Interaction of the compact ob- 
ject with the disk perturbs the disk and strips away a 
part of the disk which escapes from the system to the 
radio emission region. Gradual buildup of the equatorial 
disk due to ejection of matter from the Be star leads to 
the increase of the disk density and/or disk size. Such a 
scenario implies that extended radio emission shouldn't 
have a clear jet-like morphology, but rather have an ir- 
regular morphology varying over the orbital and superor- 
bital cycle. Such a variable morphq l ogu is indeed observed 

¥ "" 

tion of the disk could be explained by the increase 
of the time of accumulation of energy sufficient for 
ejection. The kinetic energy needed to strip away a 
part of the disk is comparable to the gravitational 
binding energy of the disk, U ~ GxM*pdRdHd 



in the radio band |Dhawan et al. ( [2006 [ ) 
The shift of the phase (fix of ejection of a por 



[R d /10 12 cmj [H d /10 12 cmj 



10 40 [pd/10 13 g • cnT 

[M*/10Mq] erg, where pd-,Pd and Hd are the density, 
radius and thickness of the equatorial disk. Such energy 
should be transmitted to the disk by the compact 
object at each disk-compact object interaction event. 
Suppose that the compact object injects energy in 
the disk at a constant rate P erg/s at each inter- 
action event. The energy sufficient for ejection of a 
part of the disk is then accumulated on a time scale 

T ej ~ U/P ~ 1 [P/10 35 erg/s]" 1 d. Accumulation of 
mass in the disk leads to the increase of T e j and, as a 
consequence, to the shift of (fix- 

If the disk size reaches the size of the binary orbit, 
the compact object always moves inside the disk and 
continuously perturbs it. Studies of the high-mass X- 
ray binaries with Be stars show that in this case the 
compact object induced instabilities in the disk might 
lead to destruction and complete loss of the disk. Our 
hypothesis is that the loss of the disk corresponds to the 
period in the superorbital phase range $ ~ 0.4—0.5 when 
the strength of X-ray and radio flares decreases and no 
systematicperiodic variabili ty of t he source is observed 
(see Fig. [Tj |3| and IGregory | ( [2002] ) ) . Ejection of matter 
from the equatorial regions of Be star leads to formation 
of a new disk at around <£ ~ 0.5 and a new cycle of disk 
growth/decay starts. 

In such a scenario, regularity of the superorbital mod- 
ulation in the system could be readily explained. The 
constant growth rate of the equatorial disk of Be star is 
determined by the stable rotation of the Be star. The 
period of superorbital modulation is determined by the 
fixed time scale on which the disk growth to the size com- 
parable to the size of the binary orbit. This scenario for 
the origin of the orbital and superorbital modulations of 
the source flux could be tested using the Ha data which 
provide a diagnostic of the state of the equatorial dis k 
of Be star flZamanov et al.||1999] |McSwain et al.||2010| ). 
Growth and decay cycles of the disk lead to the varia- 
tions of the overall strength and shape of the Ha line. 
Variability of the line intensity and profile on the time 
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sc ale of the superorbital modulation was demonstrated 
by Zamanov et al. (1999). Monitoring of the Ha line on 
the time scale ol several superorbital cycles would show 
if the the observed variability corresponds to the peri- 
odic buildup and decay of the disk. Periodic ejection of 
a part of the disk as a result of the compact object-disk 
interaction might be responsible for occurence of tran- 
sient red or blue "shoulders" of Ha line as observed by 



McSwain et al. (2010). Systematic re-observation of the 
repetition ol occurence of the shoulders in many orbital 
cycles and correlation of the phase of occurence of the 
shoulders with the phases of X-ray flares would provide 
a direct test for our model. 

Note added in proof: During t he process of p ublication 
of this article, a similar study by Li et al. (2012 ) appeared 
in press. 
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